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Abstract 

G\ ■ 

For the description of the Universe expansion, compatible with observational data, 
a model of modified gravity — Lovelock gravity with dilaton — is investigated. D- 
dimensional space with 3- and (D-4)-dimensional maximally symmetric subspaces is 
F ^ \ considered. Space without matter and space with perfect fluid are under test. In 

Lin ' various forms of the theory under way (third order without dilaton and second order 

QO ■ — Einstein-Gauss-Bonnet gravity — with dilaton and without it) stationary, power- 

CN \ law, exponential and exponent-of-exponent form cosmological solutions are obtained. 

i Last two forms include solutions which are clear to describe accelerating expansion 

^ | of 3-dimensional subspace. Also there is a set of solutions describing cosmological 

i | expansion which does not tend to isotropization in the presence of matter. 

; Introduction 

> ; 

At present time there are numerous observational data known to be incompatible with 
Standard Cosmological Model. On the one hand, accelerating expansion observations from 
supernovae type la pQ and gravitational lensing [2] allow us to calculate metric tensor. On 
the other hand, evaluating of amount of visible matter, energy-momentum tensor can be 
obtained. However, it is impossible to satisfy Einstein equations by plugging in these values. 
Then there are two possibilities: there is a great amount of invisible matter or Einstein 
J> ! equations is not true. These possibilities point out two approaches to the problem: to 
develop theories of dark matter and dark energy or to modify theory of gravity. 
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In the present article we are engaged in the second approach. At that we will give atten- 
tion exclusively to obtaining cosmological acceleration in the space with extra dimensions or 
without them, with perfect fluid or without matter, and to obtaining cosmological solutions 
which do not tend to isotropization. The latter are important for the reason that extra 
dimensions should be small and then they should not expand such as visible ones. We will 
not deal with gravitational lensing, galaxy moving in clusters, rotating curves of galaxies 
and so on in the framework of theory under investigation. Also we are not dealing with any 
issues related to quantization. 

Modified gravity has its beginning in 1920-th. The most popular theories are Brans- 
Dicke theory [3l 0], Lovelock gravity [5] and /(i?)-gravity (see, e. g. [6]). However, for a 
long time these theories were not useful for explanation of experimental data incompatible 
with general relativity. 

Here we will investigate a scalar-tensor extension of Lovelock gravity — Lovelock gravity 
with dilaton which might have its origin in low-energy limit of string theory. Lovelock 



e-mail: ikirnos@mail.ru 
2 e-mail: andre@tspu.edu.ru 



1 



gravity with dilaton contains scalar field ip (dilaton), metric tensor g^, matter fields $ 7 and 
is described (in D-dimensional space-time) by Lagrangian 

m 

c = ^^W^r^^r^™ • • •i? A2p _ lA2p <T2p - 1CT2p +^^^-^(^)+/: M (^,^), 

p=l 

(1) 

where 

m = -D, if D is even, 
2 ' 

m=-(D-l), if D is odd, 

<y p (<p), V(ip) are arbitrary functions of dilaton, 5^7.'^* is the generalized Kronecker delta 
and is equal to 1 if V\ • • • is even transposition of \i\ • • • /i^, to —1 if odd one, and to zero 
otherwise; Cm(& , 9nv) is the matter Lagrangian. We shall call terms 

r OJfAl—Aap D C1C2 p 0-30-4 T3 V2p-l<yip fc\\ 

H> — zo <ri - o'2 p - a AiA2 rl A 3 A 4 ' ' ' ""^p-^ l z J 

as Lovelock Lagrangians of p-th order. 

Theory under investigation is outstanding by the fact that its field equations are nonlinear 
with respect to second derivatives of metric tensor but does not involve higher derivatives. 

Now say a few words on the researches in the framework of the theory under consideration, 
which have already been done. Usually (see, e. g. [3 El El) only 2-nd order (i. e. a p = 
Vp > 2) of Lovelock gravity without dilaton (so-called Einstein-Gauss-Bonnet gravity) is 
investigated. Solutions for more complicated variants of the theory are not large in number. 
Investigations of third order Lovelock gravity (without dilaton) can be found in [TUl [TTJ [T2] , 
studies of the second order with dilaton can be found in [T31 [33]. Moreover, C. C. Briggs 
obtains explicit formulae for the 4-th and 5-th Lovelock tensors. Also we should draw 
attention to researches in f(R, £2)-gravity [13J, [T7] and to works in the 3-rd and 4-th orders 
of theory which is similar to Lovelock gravity and obtained from string theory low-energy 
limit as well as the latter [TSJ QIJ] . 

In the present paper a set of solutions in second order with dilaton and without it, also in 
third order without dilaton are obtained. In the most of solutions extra spatial dimensions 
are assumed to exist. Unobservability of them is explained by Kaluza-Klein approach (see, 
e. g- [20]. p. 186] and references therein) which is briefly the following: extra dimensions are 
compactified on so small scale that it is impossible to observe them (by present day device). 

In the first section seven-dimensional third order Lovelock gravity without dilaton is 
considered. Second section is devoted to the second order with dilaton in spaces with various 
number of dimensions. 

1 Third order Lovelock gravity without dilaton 

Because it seems impossible to consider general case of Lovelock gravity in a space with 
great number of dimensions, let us now discuss third order Lovelock gravity without dilaton 
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and without cosmological constant: fields are g pv and Lagrangian is 

^Loveiock3 = R + a 2 C 2 + a 3 C 3 + £ A /($ / , <7 M „). 
Here «2> a 3 are constants. General expressions for 2-nd and 3-rd Lovelock Lagrangians are: 

£ 2 = R^pR^ - ^R^ V + R 2 (3) 

is Gauss-Bonnet Lagrangian, 

f o TyuvaK p Pot 1 S E?/^ Dff« f?PT 1 O/l TDUvaK D Dp i 

i-3 — /-nr ti aK p T nr av + e/t r CT0 .ft j, t -tv UK + z4/ir K aK1/p K r ^ + 



(4) 



+ 3RR^ K R^ aK + 2AR^ vaK R oll R KV + lQR^R va R\ - YIRR^R^ + i? 3 



is third Lovelock Lagrangian. 

By variation of action S = J C-Loveiock^d D x one may obtain 



Gg! + rfg + 036$ = (5) 



where 

Gg) = ^ - (6) 

= 2(— -R M(TKT -R Kro ' v — 2R WV(J R pcJ — 2R IMT R a v +RR fiV ) — -(Rp^R^^ — 4R a/3 R a ^ + R 2 )g^ u , 

(7) 

and is written down according to [15J: 

Gp 3 } = -(^—g piJ R 3 + 12 g^RRapR 01 ^ — "&gp, v RR a p aK R a ^ aK — IQg^Rj 3 R^ R a a + 
+24gp U R a/3 R aK R aal3K + 2Ag^ y R ( J 3 R oiaKp R 0(TKp + 2g flu R a/3 r7K R aK px R 



—89uuR a 3 7K R ao aX R K \ l3p + QRuuR 2 — 2ARR °R av — 24Rr lu R Ql3 R al3 + 



+4:8R fJ a Rj 3 R/3 l/ + 48 R^R? a R a $ av + 6i? /J j / i? Q/ g (TK i? a/3crh ' — 2AR iia R v p aK R a ^ ah ' + 
+24RR paUfi R 7K + 12 RR fia /3 a R l/ al3 ' 7 — 48R fiai , l 3R (T 'R crl3 — 4&R lwl p a Rj > 'R acr + 
+48R liau/ 3R crK R ae7,3K — 24R lia p cr R Kal3cr R KV — 24R lia p a R aK R v J 3a + 
+48R fJ a P' T Rp K R aKua + 24:R fiau i3R a crKp R /3 ' 7h ' p + YIR^ 11 R^ p R Kp(XU + 



1.1 Cosmological equations 

Now consider seven-dimensionall flat space and assume metric tensor to get the form 

g, v = diag{-l, a 2 (t), a 2 (t), a 2 (t), b 2 (t), b 2 (t), b 2 (t)}. (9) 



Furthermore, let T^ v = 0. 



From such a metric one can obtain nonzero Christoffel symbols: 



T u = ad, 



1 aa — uu L iO — L Oi 



a 

j 

o 



1 cO — 1 0c 



(Latin indexes from the middle of alphabet i,j,k,... run over visible subspace, and Latin 
indexes from the beginning of alphabet a,b,c,... run over extra subspace; index notice the 
time coordinate; Greek indexes run over all the space). Nonzero components of Riemann 
tensor are 



R° 



R 



R 



tOi 



0/0 



J I.I 



a(t)a(t); 
d(t) 



R° 



cOc 



R\ 



a{ty 

d 2 (t), i^j; 
a{t). 



R 



b(t) 



a(t) 



b(t)b(t); 



0c0 " b{ty 

R c dcd = b 2 (t), c^d; 
b(t) 



(10) 



R 



b(t) 



a(t)d(t). 



Now the field equations ([5]) are 

' H 2 + 3Hh + h 2 + 12a 2 H 3 h + 3Qa 2 H 2 h 2 + 12a 2 Hh 3 - 2A0a 3 H 3 h 3 = 0, 



H (2 + 24a 2 Hh + 24a 2 h 2 - 288a 3 Hh 3 ) + h(3 + 48a 2 Hh + 12a 2 H 2 + 
+ I2a 2 h 2 - A32a 3 H 2 h 2 ) + 3H 2 + QHh + Qh 2 + 72a 2 H 2 h 2 + 
+ 72a 2 Hh 3 + 24a 2 H 3 h + 12a 2 h 4 - A32a 3 H 2 h 4 - 288a 3 H 3 h 3 = 0, (11) 

if (3 + 48a 2 Hh + 12a 2 h 2 + \2a 2 E 2 - 432a 3 H 2 h 2 ) + h{2 + 2Aa 2 Hh+ 
+ 2Aa 2 H 2 - 288a 3 H 3 h) + 3h 2 + QHh + 6H 2 + 72a 2 H 2 h 2 + 
+ 72a 2 H 3 h + 24a 2 Hh 3 + 12a 2 H 4 - A32a 3 H 4 h 2 - 288a 3 H 3 h 3 = 0. 

Here 

H(t) = a(t)/a(t), h(t) = b(t)/b(t) (12) 

are Hubble parameters for visible and extra dimensions respectively. 

Note that third equation is consequence of two other equations. A reader is asked to 
verify this correlation by himself. 

3 We elected 7-dimensional space for two reasons. First, just in such a space nonzero third order of 
Lovelock gravity arises. Second, in such a space we obtain general exact solution of the present form for the 
second order theory (see below). 
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It seems reasonable to begin consideration of system (fTTI) with allowing a 2 = = 
and solving therefore 7-dimensional Einstein equations. But all solutions of them are only 
particular cases of generalized Kasner solution (see [21] and [22} §117] for 4-dimensional 
case). Particularly, in our case (j2J) there is a solution 



a(t) 



,(o). 



1 



,(°). 



1 



(to - t) a 



3+V5 
(t - t) 6 (^+ 2 ) 

b{t) = 6(°)(* - t)T&&> « 6 (0) (to - *) a53Q 
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(13) 



where a(i) describes accelerated expansion of visible subspace. 

1.2 General solution in the second order 

Assume a 3 = in (ITT|) . Then we have 

( H 2 + 3Hh + h 2 + 12a 2 H 3 h + 3Qa 2 H 2 h 2 + \2a 2 Hh 3 = 0, 

H (2 + 2Aa 2 Hh + 2Aa 2 h 2 ) + h{3 + 48a 2 Hh + 12a 2 H 2 + 
+ \2a 2 h 2 ) + 3H 2 + 6Hh + 6h 2 + 72a 2 H 2 h 2 + 
+ 72a 2 Hh 3 + 2Aa 2 H 3 h + \2a 2 h A = 0, 



H(3 + A8a 2 Hh + 12a 2 h 2 + V2a 2 H 2 ) + h(2 + 24a 2 Hh+ 
+ 2Aa 2 H 2 ) + 3/i 2 + 6iZ7i + QH 2 + 72a 2 H 2 h 2 + 
+ 72a 2 H 3 h + 12a 2 # 4 = 0. 



(14) 



From the first equation there are 3 possibilities: 

1 „ 3 - y/h 



H 



I2a 2 h'' 



H 



h. 



H 



3 + VE 



h. 



Second and third possibilities are satisfied in all cases (H > if h < 0), first one 
a 2 > 0. Consider them one by one. 



under 



1. H 



1 



12a 2 h 



Plugging into the second equation we have 



1728a^ 6 + 1 



Then 



where x obeys equation 



12a 2 (12a 2 /i 2 + U4a 2 h 4 + 1) ' 

1 nr 



h = -W — x, 

O V « 2 



x + 3^x — 5x — 57X + 3x + 7 = 0, 
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with 



2JdT 2 



7 = tan 
3 — y/h 

2. H = — h. Plugging into the second equation we have 

• _ 3h 2 (-A0a 2 h 2 + 16a 2 h 2 V5 - \fh + 5) 
2(-y/E + l8a 2 h 2 V$ - 30a 2 h 2 ) 

Then h obeys equation 

8 arctanh (—yfa^h + \/5a 2 h) yfa^h + 1 — 6(t + t )h + \fh = 0. 



3 — \~ \/~^) 

3. H = h. Plugging into the second equation we have 

2 

■ h _ 3h 2 {tf$u 2 h 2 - 5 - V5 + 16a 2 h 2 V5) 
2(30a 2 h 2 + l8a 2 h 2 V5 - V$) 

Then h obeys equation 

8 arctanh U/a^h + y/ba^h) yfa^h — 1 + 6(t + t )h + v5 = 0. 



In all cases above there is only parametric dependence H(t) and h(t). Some explicit 
solutions in Einstein-Gauss-Bonnet gravity for 7 or other dimensions will be obtained in the 
second section. 

1.3 Exponential solution in the third order 

Now consider equations of the third order Lovelock gravity with constant Hubble parameters: 

H = 0, h = 0. 
Then the system (fTTl) takes a form 

' H 2 + 3Hh + h 2 + 12a 2 H 3 h + 36a 2 H 2 h 2 + 12a 2 Hh 3 - 2A0a 3 H 3 h 3 = 0, 

3H 2 + 6Hh + Qh 2 + 72a 2 H 2 h 2 + 72a 2 Hh 3 + 2Aa 2 H 3 h + \2a 2 tf- 

- 432a 3 # 2 /i 4 - 288a 3 H 3 h 3 = 0, (15) 

3h 2 + QHh + QH 2 + 72a 2 H 2 h 2 + 72a 2 H 3 h + 24a 2 Hh 3 + 12a 2 H 4 - 

- 432a 3 H 4 h 2 - 288a 3 H 3 h 3 = 0. 

Subtracting the second equation from the third one we have 

(H 2 - h 2 ) [l + \Qa 2 Hh + Aa 2 (H 2 + h 2 ) - UAa 3 H 2 h 2 ] = 0. 
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Then h = —H is a solution of this equation. 

Plugging this equality into the second equation of (|T5l) one can obtain 4 various solutions: 



, —a? ± \l «n + 12a 3 
H = ±\ V 2 -. 16 

Because of the expansion of visible dimensions let us take only H > 0. Then we have 



, -q 2 ± yjaj + 12q 3 / -q 2 ± ^/oj + 12q 3 

~ V/ 24a 3 ' V 24 "3 ' 



Plugging /i = —H into the first equation of (1151) we obtain 



„ = ± i -a 2 ± yg + (20/3)a 3 (lg) 



Comparing that with (fl6|) we have 



Then 



40a 3 



1 .2 



«3 = ~Y2 a 2- 



V2tt2 V 2 «2 

Therefore (because of (TP2]) ). 



a(t) = c ' exp (vS;J' 6(t) = C2exp l-7S;J' (19) 

where Ci, C2 are arbitrary positive constants. 

It is clear that a(t) > 0, a(t) > 0, i. e. the abovementioned solution describes accelerated 
expansion of visible dimensions. At that extra dimensions shrink. Then it is possible that 
visible and extra dimensions were equivalent and the Universe look as 4-dimensional one 
only after expansion of one subspace and contraction of another one. Solution ffl9|) seems to 
be useful for the description of inflation. 

Finally, if T^ v = 0, a 2 > and a 3 = — ^7^2 then system (JSJ) has solution (j5J) where 
functions a(t) and b(t) are expressed by ( 1T91) with arbitrary positive constants C\ and C 2 - 



2 Einstein-Gauss-Bonnet gravity with dilaton 

Such a theory is the following: fields g^ u , <p, Lagrangian 

C EGBd = R + g^dtfdw - V(<p) + e(cp)C 2 + £ M ($ 7 , g^). (20) 

Here s((p), V((p) are functions of dilaton <p, C 2 is a Gauss-Bonnet Lagrangian (j3J). Theory 
under consideration is different from generalized Brans-Dicke theory (see, e. g. [2U]) even 
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in 4-dimensional space, that is why we investigate both (3+l)-dimensional space and spaces 
with extra dimensions (where Einstein- Gauss-Bonnet gravity without dilaton is sensible). 
Variating the action with Lagrangian (J20l) we get field equations: 

'11 1 

R^u - -Rg^u - -g»ug a(3 d a (pdf3<p + -gp,V((p) + e((p)G® + d^tpd^- 

- 4 - l ~Rg^ Ueiip) - 4 {r^ v + R^g^ V„V /J e(^)+ 



+ 8i2%V Q V„e(p) - 2RV ll V v e(<P) 



2L> + V'(<p) - e'(<p)£ 2 = 0. 



8ttG 



(21) 



T 



„4 A*"' 



Here is the second Lovelock tensor ([7]). 



2.1 Cosmological equations 

Consider space of/} = p + q + 1 dimensions with two maximally symmetric subspaces: 
j9-dimensional and g-dimensional. Square interval in such a space is 



ds > = -e 2 Mt) dt 2 + e 2m(t) ds 2 + e 2u 2 (t) ds 2 



(22) 



where ds 2 and (is^ are square intervals in p- and g-dimensional subspaces respectively, 
u (t) , Ui(t) , u 2 (t) are arbitrary functions of time t. 

If metric is (1221) then Christoffel symbols are (as above, Latin indexes from the middle 
of alphabet i,j,k,... run over visible p-subspace, and Latin indexes from the beginning of 
alphabet a,b,c,... run over extra g-subspace; index notice the time coordinate; Greek 
indexes run over all the space) 



r° 



oo — M o, 



T° ab = u 2 e- 2uo g ab , 



pi pi 

1 jk — 1 jkt 



pa pa 

1 be ~ 1 bci 



rna — r"n — v>2. 



(23) 



Riemann tensor, Ricci tensor and scalar curvature are 



R°io 



-2uq 



Xg 



R° 



aOb 



-2uo 



Yg 



R l jki = e ~ 2u °A p (5 l k gji — Slgjk), R t a jb 
R a bcd = e- 2u «A q {5 a c9bd - S a d g bc ), R 00 



ab ■ 



e 2ua uiu 2 5)g abl 
-pX - qY, 



R i3 = e~ 2u ° {X + {p-l)Ap + qinu 2 ) g ih R, 



ab 



-2u 



(Y + (q- l)A q + puiu 2 ) g a b, 



R = e~ 2uo [2pX + 2qY + Pl A p + q x A q + 2pqu 1 u 2 ] . 
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Gauss-Bonnet Lagrangian (j3J): 



C 2 = e- 4 " \p 3 A 2 p + 2 Piqi A p A q + q 3 A 2 q + Au 1 u 2 (p 2 qA p + pq 2 A q )+ Ap iqi u 2 u 2 + 

+ 4pX[(p - l) 2 Ap + qi A q + 2(p - l)guiu 2 ] + 4gl>iA> + (? - 1)2^ + Ml ~ 1)^2] j 



Here we introduce the following notations: 

A p = iij + v P e 2iu °- Ul) , A q = ii 2 + a q e 2 ^ u °- U2 \ 
X = u x — uqUx +u(, Y = u 2 — UqU 2 + u% 

(p — m) n = (p — m)(p — m — l)(p — m — 2) . . . (p — n), 
(q — m) n = (q — m)(q — m — l)(q — m — 2) . . . (q — n), 



(24) 



(25) 



a ^W^) a ^W^ry (26) 

where R p R q are internal curvatures of p- and q- dimensional subspaces respectively, Y % - k and 

r^ c are internal Christoffel symbols. 
Now field equations are 



- Pl A p + X -q x A q +pqii x u 2 + hp 2 - l -e- 2u »V{y) + h~ 2u °e( V ){p 3 A 2 p + 2p iqi A p A q + 

(27) 

+ q 3 A\ + Au x ii 2 (p 2 qA p + pq 2 A q ) + 4p iqi u 2 ii 2 } + 2e- 2u °e'{^{A p {p 2 qu 1 + y ' 

+ qp\u 2 ) + A^pqxUx + q 2 ii 2 ) + 2u 1 u 2 (piqui + pqiu 2 )} = ^f-T 00 , 
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e- 2 "»{(l - p )X-qY--(p- l) 2 A p - -q x A q - (p - l)g«i« 2 }^- + -e^^g^ 
+ \v{y)gij - l -s{y)e-* u ° gij {(p - l) 4 A 2 p + 4(p - l) 2 gi«?«i + 4(p - 1) 3 A P X+ 
+ 4(p - l) 3 gMiM 2 A p + 8(p - l)gitii«2y + 4(p - l)g 2 «iW2^ + 4(p - l)giA,X+ 
+ 8(p - l) 2 gMiM 2 X + 4(p - l) 2 gA p F + 2(p - l) 2qi A p A q + + 4g 2 A g F}+ 
+ 2e- 4M( %{-(eV + e'0 - £ , u <p)[(p - l) 2 A p + q x A q + 2(p - l)g«iu 2 ] + 
+ e'« 1 ^[-2(p - 1) 2 X - 2(p - l)gY - (p - 1) 3 A P - (p - \)q x A q - 

- 2(p - l) 2 guiu 2 ] + £'« 2 0[-2(p - l)gX - 2 qi Y - (p - l) 2 qA p - 

i-. 87rG ^ 

- g 2 A, - 2(p - l)giwiw 2 j} = 

e- 2 ««{(l - q)Y - pX - -{q - l) 2 A q - X -p x A p - (q - ^pii^hab + ^e~ 2u ^ 2 g ab + 
+ l -V{ V )g ab - \e( V )e- 4uo g ab {(q - l) 4 A 2 q + 4(g - \) 2 p x u\u\ + 4(g - 1) 3 ^Y + 
+ 4(g - l) 3 pw 1 w 2 A 9 + 8(g - l)pi«iW 2 X + 4(g - l)p 2 wiw 2 ,4 p + 4(g - l)piA p F+ 
+ 8(q - l) 2 puiU 2 Y + 4(g - l) 2 pA,X + 2(g - l) 2Pl A p A q + p 3 ^ + 4p 2 A p X}+ 
+ 2e- 4uo ^ a6 {-(eV + e'0 - e'u cp)[(q - l) 2 A g + Pl A p + 2(q - l)pu 1 u 2 ]+ 
+ £ '<w[-2(g - l) 2 y - 2(g - l)pX - (g - 1) 3 A 9 - (g - 

- 2(g - 1)2^^2] + e'u 1 ^[-2(g - l)pY - 2 Pl X - (g - l) 2 pA, - 

- p 2 A p - 2(g - l)pi«i« 2 ]} = -^r T ab, 
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2\(p + (-mo + pin + qu 2 )ip} - V'{(p) + e\ip)e~ 2uo ^p 3 A 2 p + 2p 1 q 1 A p A q + q 3 A 2 q + 

+ Au 1 u 2 (p 2 qA p + pq 2 A q ) + Ap iqi u\ul + ApX[(p - l) 2 A p + qi A q + (30) 

+ 2(p - l)guiU2] + 4gl> 1 A p + (g - l) 2 A g + 2p(g - l) Wl w 2 ] 1 = 0. 

These equations are equivalent to those in [T4| if we substitute g^d^tpdytp by —-g^d^ipdvip 
in Lagrangian and put e((p) = a 2 e~ lip , V((p) = 0, = 0. 

Henceforth, we may put uq = (for simplification), p = 3 (to identify p-subspace with 
visible space). 

2.2 Stationary solutions 

Let us now turn to find solutions of (|2T1) - (I3"0"I) under p = 3, Mo = 0. The simplest solutions 
are stationary ones. Hence put 

Ui = const, u 2 = const, <p = const . 

Then ([24]) are 

A p = a p e~ 2ui , A q = a q e~ 2u2 , X = Y = 0. 

Consider space with homogeneous dust, i. e. T o ^ (other T^ u = 0). 
After that, system (I2TT ) ([30]) get the form of algebraic equations 



3A p + X -q x A q - ±V(<p) + h{ip){l2 qi A p A q + q 3 A 2 q } = ^T 00 ; (31) 



-A p - l - qi A q + ±V{<p) - l -e{ip){A qi A q A p + q 3 A 2 q } = 0; (32) 



-(q - l) 2 A q - 3A P + -V(<p) - -e(ip){(q - l\A 2 q + \2{q - l) 2 A p A q } = 0; (33) 



V'(ip) - e'(<p){l2 qi A p A q + q 3 A 2 } = 0. (34) 
From fl34l) and (l3Tj) we have 

34, + \qiA q = \V{ V ) - + ^oo- (35) 
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Try to find linear combination of (|3T]) - (!33|) in order to cancel terms with e(ip). Let a, /3, 7 
are coefficients for (|32|) . (1331) and (}3~Tj) in such a combination. Then we need 

aq 3 + (3(q - 1) 4 = 7<?3, 
4^ + 12/5(^-1)2 = 127^. 

Therefore 

7 = ( 1 - -1/3, a = -p. 



qj q 

Now put P = q. Then a = 3, 7 = g — 1. 

Now multiplying (1321 by 3/ (g + 1), ( 1331) by g/ (g + 1) and putting them together we have 
(taking (1341) into account) 

From this and fl35|) one can get 

>) + 4t^)^ - ^ - 0. ,37) 
g + 1 g + 1 £'(</?) c 4 

Put now 

V(tp) = ae~ av , e(<p) = be~^. 

Then (f3~Tj) get the form 

1 / cA „,„ 87rG? 

1 + - ] ae~ av -T 00 = 0. 



q + l\ ' PJ ~ c 4 

It is easy to see: 



1 f P(q+l) 8nG 

<P = ln i / , m • -^ T oo t • 38 

a {a(a + P) c 4 



Plugging this if into (1351) we have: 



Plugging A„ derived from that into ( 1341) we can get 



« - 2)(, - 3) - 2,(, - 1)} + + = « I «£±1) 1 

y a + p po [ a(a + p) J 



a-f3 



where 

_ SvrG 

^ = — T--MX)- 

C 4 

That is quadratic equation on A q , which solutions are 



A - = 2^W) • < 4 °> 

12 



where 

2 



Then, taking 0391) into account, we obtain 



(41) 



Finally, 



A p = lx (1 " 9)a + ( J + 3)/3 - (42) 
p 6 a + /3 6 q K J 

u x = -\\n^, M2 = -iln^ (43) 
l o v l a q 

and (p is ( 1551) . 

It should be noted that the following constraints was applied: 

Tqo^O, a^O, b^O, a ± 0, /3 ^ 0, 

a^-/9, g^O, g^2, q ^ 4, a p ^ 0, a, ^ 0. 

It is easy to derive solution for (3+l)-dimensional space with perfect fluid of arbitrary 
equation of state parameter. At that we should not specify V(ip) and e((p) because of s(ip) do 
not participate in equations and V((p) is specified from those. Dilaton also do not contribute 
in equations, therefore we should solve just Einstein equations with cosmological constant. 
Solution is 

87rG_ 1 fl+w8nG r 



V = (1 + 3w)^ r T 00 , m = -- In — —T 00 

C I \ lOp c 

Here u> is equation of state parameter (p = we, p is pressure, e = T 00 is energy density). A 
particular case (when w = 0) was derived by A. Einstein in 1917 [23] . 



2.3 Exponential solutions 

For the dynamical solutions we need to do further simplification of (T2"T|) - (l3"0l) . Therefore, in 
addition to p = 3 and u = 0, put a p = a q = i. e. subspaces are flat. In addition to 
simplicity, such a condition is caused by Cosmic Microwave Background observations [24"1 [25] 
indicate the flatness of visible subspace. For extra subspace a q = is only a simplification. 
After that, equations (I27 1) -( l30l) are 



3 «l + + 3gitiM 2 + - -V((p) + ^£(v3){36gi?iiM2 + 93«2 + V2u x ix 2 {2qu\ + q 2 u 2 2 )} + 

+ 2e'(p{Qqu 3 1 + l%qu\u 2 + 9q 1 u 1 ii 2 l + q 2 ul) = ^^T 00 . 

(44) 
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1 11 

{-2«i - qu 2 - Suj - -q(q + l)u\ - 2qu 1 ii 2 }g ij + -0 2 ftj + -V^y?)^-- 



(45) 



- ^£(v?)5ij{8iiiu 2 (giw 2 + 2qui) + 4w 2 (4giWiw 2 + 2qu\ + g^) + (? + 1)<?2« 2 + 
+ 4g(5g - 2>)u\u 2 2 + 8ggiUiui + 16gw?w 2 } + 2g ij {-(e"ip 2 + e'(p)[2u\ + q x u\ + 



+ 4g-ui-u 2 ] + E'u\(p\-A.U\ - 4qu 2 - 4u\ - 2(g + l)qu 2 2 - Aqu l u 2 } + 

8tvG 

+ e'ii 2 ^\-4qu x - 2q x u 2 - 6qii 2 - qq x ii\ - 4g 1 w 1 w 2 ]} = ~^ T ij- 



{-3ui - (q - l)u 2 - 6ul - ^q\u\ - 3(g - l)uiu 2 }g a b + \^9ab + ^V(<f)g ab - 

- '^e{ip)g ab {12u 1 {2ij^ + (q- l) 2 u\ + 4(g - l)«itt 2 ) + 4ii 2 (6(g - l) 2 Mi« 2 + 
+ (g - 1) 3 «| + 6(g - l)u?) + 24w^ + g 3 w| + 24(g - l)(2q - 3)u?u|+ 
+ 72(g - l)uf« 2 + 12(g - l)(g - l) 2 «i^} + 2^ a6 {-(eV + e'$[6u?+ 
+ (g - l) 2 ul + 6(g - l)«iU2] + £'w 2 y?[-6(g - l)«i - 2(g - l) 2 ii 2 - 

- 12(g - l)u? - (q - l)(q - \) 2 ii\ - 6(g - l) 2 «i« 2 ] + e'ui0[-12iii- 

- 6(g - l)u 2 - 18«? - 3^ - 12(g - l)«i« 2 ]} = ^T ab . 

-2</3 - 2(3tti + qu 2 )(p + V"(<p) - e'(<p){12wi[2u? + g^ + 4gM!M 2 ] + 4gii 2 [6M? + (g - 1) 2 «2+ 

+ 6(g - l)uiU2] + 24iif + (q + l)g 2 i^ + 24g(2g - l)u?u| + 72qii\u 2 + 12ggiUi«|} = o. 

(47) 

Find at first solutions without dilaton, without matter and with constant Hubble param- 
eters: 

(p = 0, V((p) = 0, T^ u = 0, iii = const, ii 2 = const . (48) 



14 



(46) 



Then system (jUD-flHH) will be a system of algebraic equations for which we have found two 
analytical solutions for arbitrary q and negative e: 



ui = U2 = ± ■ (49) 

-q(q+l)e 



Also particular cases from q = 1 to q = 22 (22+4=26 is required for the bosonic strings, 
and in the case of q = Lovelock gravity is just Einstein gravity) have been studied, but 
solutions different from (|4"9~]1 have been obtained only for q — 3 (i. e. for 7-dimensional space 
just as in section 1): 



, , 3 + Ci-Jl 575 + 257(2 *Jl . , , 3 + 

"■ = -61/ - 3010+ 1346^ ' " 2 = Cl V^^' <50) 



and 



,/ 3 + C 2 V5 . ,/ 3 + C2V / 5 575 + 257C 2 \/5 

= -"-f'V^ - ' 3010 + 13466^ <51) 

where constants Ci an d C2 take values of +1 and —1 independently from each other, and 
e > 0. 

Therefore scale factors are 

a(t) = e ui = aoe" 14 , b(t) = e U2 = b^ 1 . 

It is clear that solutions ( )49l) are not useful for us by the following cause: when visible 
subspace expands, extra subspace expands too, then extra subspace must be visible in this 
case. But solutions ( 1501) for (1 = — 1 and ( 15T1) for (1 = +1 satisfy our purpose. 

Now let us try to obtain exponential solutions in the presence of perfect fluid. For that 
substitute conditions (j4"8l) by 

(p = 0, V(ip) = 0, Too = e, Tij = wegij, T ab = weg ab , iii = const, u 2 = const . 

(52) 

After plugging those into fj44l) - (|46j) and subtracting factors and g ab we see that left-hand 
sides of equations are independent of time. Hence the right-hand sides also must be constant. 

From 00-component of local conservation law for energy-momentum tensor (V M T M o = 0) 
one can obtain (taking (1521) into account) 

e = e exp[-(l + w)(3u 1 + qu 2 )t\. 

Therefore e = const under at least one of a two conditions (here H = iii, h = U2)' 

1. w= -1; 

2. h = --H. 

Q 

In the first case matter can be described by cosmological constant, in the second one 
comoving bulk is constant. In the latter case equations (S3D — (SSD , as equations on H,e,w, 
have two solutions: 
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1. 

H = 0, e = 0, w is arbitrary, 
i. e. flat space with Lorenz metric. 

2. 

3c 4 H 2 (-q 3 - 3q 2 + 3eH 2 q 3 + 54eH 2 q 2 + 81eH 2 q - 162eH 2 ) 



H is arbitrary, 



w 



16vrGg 3 

eH 2 q 2 -q 2 + 15eH 2 q - 18eH 2 



-q 2 + 3eH 2 q 2 + A5eH 2 q - 54eH 2 ' 

(53) 

i. e. one can obtain any value for H by matching energy density e and EoS parameter 
w. It is clear that h < if H > 0, that's why this solution satisfies all requirements. 
Finally, such a solution describes anisotropic expansion of the Universe with matter 
which not tends to isotropization. In Einstein gravity it is possible only for maximally 
stiff fluid: w — 1. 

Now turn to the cosmological constant case: w = —1. Then it is possible to consider 
as equations on H, h and e. These have the following solutions: 

1. 

H = h and is arbitrary, 



c 4 h 2 (Q + q 2 + 5q + eq 4 h 2 + Qeq 3 h 2 + \\h 2 eq 2 + Qh 2 eq) 

16ttG ' 



(54) 



.is arbitrary, ^-( g -l)B«^ 



c 4 



2s 

(-96He 2 q 2 h 3 + 48He 2 qh 3 - 24Heqh + 48He 2 q 3 h 3 + 24Hhe + 2h 2 eq 2 + 
+ 10h 2 eq + Ue 2 q 4 h 4 - Q0e 2 q 3 h i + 82e 2 q 2 h A - 36e 2 g/i 4 - 3 - 12h 2 e) . 

(55) 

It is evident that the first solution is unsatisfactory. The second one is adequate under 
h < 0, H > 0. Such conditions are fulfilled in three cases: 

1. e > 0, h < — , sign in expression for H is arbitrary. 

^eq{q - 1) 

2. e < 0, sign is 

3. e < 0, /i < 1 ,signis"+". 

V-2e(g-2)(g-l)' 

Here we also should emphasize an existence of solutions with matter which do not tend 
to isotropization. 
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2.4 Exponent-of-exponent form solutions 

Try to obtain solutions with a dynamical dilaton. For that purpose consider equations ( 1441) - 
()47l) and notice that functions Ui(t) and v,2(t) make contribution only through derivatives 
iii, u 2 , ill, u 2 , but ip(t) participate explicitly. To find solutions with constant derivatives 
let's eliminate (p(t) by introducing new time variable. At first we put 

e(<p)=Pe-™, V^)=ae^, T, u = 0. (56) 

Now turn from time t to new variable r: 

d T = e~^ /2 d t . 

Derivatives with respect to r will be denoted by the prime '. After such a substitution and 
putting u" = «2 = V 9 " = one can obtain 

q/2,1 /2 |Q / / . 1 /2 1 . l fl r oc l2 ,2 , ,4. 

3«i + 2^ lM 2 + 3gw 1 w 2 + -v? - -a + -PiZGqiU! u 2 + q 3 u 2 + 

(57) 

+ \2v! x v! 2 {2qu! 2 + q 2 u' 2 2 )} - 2(3^'{18qu[ 2 u' 2 + % x vl x %}! 2 + 6qu[ 3 + q 2 u' 2 } = 0; 



- 3m 'i 2 - ^(5 + l )q u 2 2 - IV (u'i + - 2g«iM2 + + f - ^{-4?(3g + l)7^'wiw 2 2 - 
- 28g7<A / 1 2 ?4 - 2(g + 2)g l7 v?V 2 3 - 47 VV 2 - 2g l7 V' 2 m 2 2 + 4 <?(5g ~ 3)u[ 2 u' 2 2 + 
+ %qqiu' x u' 2 + 16qu[ 3 u' 2 + (q + l)g 2 w' 2 4 - lfryy/ui 3 - 8q-f 2 (p' 2 u[u' 2 } = 0; 

(58) 

1 ,2 R ,2 q — 1 // 3 , , , , 1 /2 a 

--gi« 2 -6% — 7^ m 2 - -7^ % - 3(g - l)u x u 2 + -<p +7^- 

- \(3{-2{q + l)(q - 1) 27 <A 2 3 - 6(g - l)(3g - 2) 1V 'u' x u' 2 - 60(g - 1) 7 ^V 1 V 2 - 

- 607^'m; 3 - 2(q - l) 2 7V 2 w 2 2 ~ V2rftpt\'* + g 3 « 2 4 + 24(g - l)(2g - 3)u' x 2 u' 2 2 + 
+ 12(q - l)(q - l) 2 u'iu' 2 + 72(g - l)u'*u' 2 + 2Au' x - I2(q - l)j 2 ip' 2 u' x u' 2 } = 0; 

(59) 

- 7 y/ 2 - 2(3*4 + qu' 2 )(f' + (3^{12q(4q - 2)w' 1 V 2 2 + (q + l)q 2 u' 2 A + 12qu' x v! 2 + 12gg 1 w' 1 w / 2 3 + 
+ 18gi7^'M>' 2 2 + 36g7^V 1 2 M , 2 + 12^'u' x + 24m; 4 + 2q 2 ^'u' 2 } + 07 = 0. 

(60) 
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Now assume that we have obtained some quantities u^u^f' which satisfy these equa- 
tions. What should they be to describe accelerating expansion of visible subspace and 
contraction of extra one? It is easy to see that the scale factor of visible subspace would be 

a(t) = e" 1 = a exp j^e^/ 2 } , (61) 

where ao, Co are arbitrary positive constants. Then its first and second derivatives with 
respect to time t would be 

a(t) = oou'jcoe^exp (^ e ^ 2 



f>> 

a(t) = -aou'^'e^exp Z^e^'M + a ^^t exp 



It is clear that both derivatives would be positive (as need for accelerating expansion) if 
u[ > 0, ff > (the latter is not necessary but is sufficient). 

By the same manner we obtain the scale factor for extra dimensions: 

b(t) = e u * = 6 exp j^ e ^ 2 j (62) 

(bo is positive constant) and its first derivative: 

b(t) = b Q u' 2 c e^ 2 exp l^e^/z 

which would be negative (as need for contraction) under u' 2 < 0. 

Therefore it is necessary to find solutions of (|57j) -( l60l) satisfied conditions 

u[ > 0, u 2 < 0, ff > 0. 

Numerical calculations give us solutions for different dimensions from q — 1 to q — 20. For 
example, 

q = 1 ; a = 1, /3 = 1, 7 = 1, f' = 0.383, u[ = 0.378, u' 2 = -1.32. 

q = 2, a = 0.01, (3 = 100, 7 = 0.01, iff = 5.09 ■ 10" 5 , u[ = 0.0175, u' 2 = -0.0803. 
q = 3, a = 0.001, (3 = 100, 7 = 0.001, f' = 2.33 ■ 10" 4 , u' x = 0.0306, w 2 = -0.0788. 

Finally, field equations in the case of flat subspaces without matter and with V(ip) and 
s((f) in the form of ( 1561) have exponent-of-exponent form solutions ( 16T1) . ( |62i) with above- 
mentioned parameters. 

2.5 Power-law solutions 

Now consider space with dust-like matter: 

T 00 7^ 0, other = 0. 
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And try to obtain solutions of system fH4l) - (|471) with scale factors of power-law form: 

a(t) = = (J-J , b{t) = = (j-Y . (63) 

Then all Einstein terms in left-hand sides of equations (|44j) - fH6l) are proportional to 1/t 2 
and Gauss-Bonnet terms are proportional to 1/t 4 . Hence solutions with scale factors (1631) 
are possible if 

Too oc — , (64) 

2 oc -, (65) 
e{ip) oc t 2 , (66) 
V(<p) oc i (67) 

Consider these conditions one by one. 

It is possible to derive time dependence of energy density from conservation law of energy- 
momentum tensor: 

T 00 = const .t' 3n - qm . 
From comparison this expression with (IM)) we have 

2 - 3n 



m 



q 



Note that under this condition extra subspace contracts (m < 0) if visible subspace ex- 
panses accelerative {n > 1) (but we haven't obtain such a solution, see below). For (3+1)- 
dimensional space condition ( |64l) leads to n = 2/3 i. e. to Friedmann solution. 

From condition ( 1651) it is easy to obtain ip{t) = -^ln(t/t 3 ), where i/;, t 3 are arbitrary 
constants (t 3 > 0). In order to avoid unnecessary complication put ip = I. Therefore 

ipit) = In f 1 



From comparison that with ( |66l) . (|67l) we see: 

V(y») = ae~ 2 ^ e(<p) = /3e 2 ^, 

where a, f3 are constants. 

Plugging all those into equations (j4"l"]) - (P4"T|) and putting 

a — > a = at 2 , /3 — > /3 = /5/t 2 , 

one can obtain a system of algebraic equations on n with parameters a, /3, g and x = 
— —T 00 (t )tl, where T 00 (t ) is energy density at some time moment t . Such a system of 
equations has solutions not at all values of parameters. Considering dimensions from q = 1 
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to q = 22 we have obtained a set of solutions with arbitrary x. Here a and (3 are functions 
of x and < n = m < 1 i. e. visible and extra subspaces expands with deceleration (and 
with the same velocity). In another set of solutions x possesses fixed values. Here in two 
cases n = m and in another cases n ^ m. These are solutions for q — 1, 2, 6, 9, 12: 



x is arbitrary, 



3 29 ,25 

a = -x , p = — x , 

7 14' M 21 21' 



n 



m = — 



x 



53 
40' 



a 



973 
'640' 







25 
'272' 



n 



3 

To' 



m 



11 
20' 



g = 6, x is arbitrary, a 



112 5077 
253 X ~2277' 



P 



729 
16192' 



-x- 



2025 
16192' 



2 
9" 



9 = 6, 



x 



51397 
39528' 



a 



21035 
13176' 







243 
3904' 



n 



27' 



5 

27' 



g = 9, 



9 



12, 



x 



x 



17 

y 

43 
15' 



a 



a 



-1, 



-1, 



= 0, 



= 0, 



n 



1 

6' 

2 

15' 



Note that in the last two cases = 0, therefore these solutions are solutions in Brans- 
Dicke theory (i. e. theory with Lagrangian (I2"0"j) without Gauss-Bonnet term). 

In all obtained power-law solutions 0<n<l,0<m<l, that's why such solutions don't 
describe accelerated expansion of visible space or contraction of extra dimensions. However 
solutions with n ^ m are interesting for another cause. In Einstein theory there is no 
anisotropic power-law solution in the presence of dust. However, in Einstein-Gauss-Bonnet 
theory with dilaton that is possible. 

Conclusion 

Different variants of Lovelock gravity with dilaton were considered in D- dimensional space 
with two maximally symmetric subspaces: 3-dimensional and (D — 4)-dimensional. Absence 
of matter and existence of perfect fluid were investigated. We have several types of obtained 
solutions: 



1. Stationary. 

2. Power- law. 

3. Exponential. 
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4. Exponent-of-exponent form solutions. 

Among the last two forms solutions which describe accelerating expansion of 3- dimensional 
subspace and contraction of (D — 4)-dimensional one were elected. Unobservability of the 
latter subspace was justified on the basis of Kaluza-Klein approach. Also a set of anisotropic 
solutions which do not tend to isotropization in the presence of matter, in contrast to Einstein 
gravity, have been obtained. Such a possibility is of importance because it allows us to assume 
that extra dimensions become small during the Universe evolution. This issue we are going 
to investigate in more detail in another work. Moreover, it would be interesting to extend the 
results of this work for account of third-order Lovelock terms. This will be done elsewhere. 

Studying of future singularities in such models would also be important. For 4-dimensional 
modified gravities this problem was considered in [261 127] . 

Unfortunately most of solutions describe only flat maximally symmetric subspaces. For 
curved subspaces there are only stationary solutions. Those are of interest as exact solutions 
of very complicated equations and as possible basis for numerical dynamical solutions in the 
case of curved subspaces. 
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